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Using semiclassical WKB-methods, we calculate the pair-production rate in the presence of 
two external fields, a strong (space- or time-dependent) classical field and a monochromatic 
electromagnetic wave. Wc discuss Pauli-supprcssion and Bose-cnhancement of the rate in 
Ql^ the presence of thermal electrons and bosons at a given temperature and chemical potential. 

Using our rate formula, we calculate the rate enhancement due to laser and thermal photons, 
and discuss the possibility that the enhancement may be significant in a plasma of electrons 
O | and protons with self-focusing properties. 
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I. INTRODUCTION 

The creation of electron-positron pairs from the vacuum by an external uniform electric field in 

n n 

spacetime was first studied by Sauter [1[ as a quantum tunneling process. Heisenberg and Euler [2[ 
extended his result by calculating an effective Lagrangian from the Dirac theory for electrons in 
a constant electromagnetic field. A more elegant reformulation was given by Schwinger based 
on Quantum Electrodynamics (QED), where the result is obtained from a one-loop calculation of 
the electron field in a constant electromagentic field yielding an effective action. A detailed review 
and relevant references can be found in Refs. 4] and 

The rate of pair-production may be split into an exponential and a pre-exponential factor. The 
exponent is determined by the classical trajectory of the tunneling particle in imaginary time which 
has the smallest action. It plays the same role as the activation energy in a Boltzmann factor with 
a "temperature" fi. The pre-exponential factor is determined by the quantum fluctuations of the 
path around that trajectory. At the semiclassical level, the latter is obtained from the functional 



determinant of t 



classical paths jfj]. An efficient technique for doing this is based on the WKB wave functions, 
another on solving the Heisenberg equations of motion for the position operator in the external 



re quadratic fluctuations. It can be calculated in closed form only for a few 



field Q]. If the electric field depends only on time, both exponential and pre-exponential factors 
were approximately computed by Brezin and Itzykson by applying Schwinger's method to a purely 
periodic field E(t) = Eocoscaot [?]. The result was generalized by Popov with a First-quantized 
calculation in Ref. [8(] to a general time-dependent fields E(t). An alternative approach to the same 
problems was more recently employed using the so-called worldline formalism [9|] , sometimes called 
the "string-inspired formalism". This formalism is closely related to Feynman's orbital view of 
the propagators of quantum fields. The functional determinant of the electron field in Schwinger's 
approach is calculated as a relativistic path integral over all fluctuating orbits of an electron in the 
external field as described in the textbook (|. In the path integral formalism the tunneling problem 
has a standard formulation and the pre-exponential factor is calculated via an orbital fluctuation 



determinant for whose calculation simple formulas have been developed in Ref. [Ql. These formula 

n n 

were evaluated by Dunne and Schubert [10|] and Dunne et al. [Ill ] for various field configurations, 
such as the single-pulse field with a temporal Sauter shape oc 1/ cosh 2uit. 

In our previous article [ijj], we have derived a general expression for the pair-production rate in 
nonuniform electric fields E(z) pointing in the z-direction and varying only along this direction. A 
simple variable change in all formulas has led to results for electric fields depending on time rather 
than space. 

The relevant critical field where copious pair production sets in is E c = m2/e = 1.3 x 10 16 V/cm, 
field intensity I c = E c 2 ~ 4.3 x 10 29 VF/cm2. For electric fields E <C E c , the pair-production rate 
is exponentially reduced by a factor exp — ttE c /E. In the laboratory, the electric field intensity 
I c is, unfortunately, extremely difficult to reach in the laboratory [lj, Q]- Motivated by these 
difficulties, people have studied possibilities of a dynamical enhancement of pair-production rate 
by time-dependent oscillating or pulse electric fields [15I. [r]]. 

One possibility is to consider the superposition of a strong but slow field pulse and a weak 



but fast field pulse, which can lead to a significant enhancement of pair-production rate [15|, ]lj 
Another is a catalysis mechanism of pair productions that has been studied in Ref. [lfj]. The setup 
is a superposition of a plane-wave X-ray probe beam with a strongly focused optical laser pulse. 
Namely, the optical laser pulse beams are focused onto a spot to yield a strong stationary electric 
field E, and the X-ray laser propagates through the focusing spot of optical laser beams. Since 
X-ray laser wavelength (frequency) is much smaller (larger) than the optical one, namely the size of 
the focusing spot, the electric field created by focusing two optical laser beams can be approximated 
by a constant classical electric field in space and time. In that spot, a large number of coherent 
photons (X-ray laser) collides with virtual pairs of the vacuum in a strong classical electric field 
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(optical intense pulse), and in consequence the pair-production rate must be enhanced. 



121 ] , by calculating the 



In this article, we continue and extend our semiclassical WKB-approach 
enhanced pair-production rate in the superposition of a strong (space- or time-dependent) classical 
field and an electromagnetic plane wave. In Sections [II] and IIIIl we present a general expression 
for the rate with a general enhancement factor. In Section HVl we apply this general expression to 
two cases; (i) a constant electric field in finite spatial region which drops sharply to zero at the 
boundary; (ii) a softened version of this, where the production takes place in a Sauter electric step 
field. In Section [V] we extend our general formalism by calculating the enhancement factor in the 
presence of coherent laser photons and thermal photons at a finite temperature. In Sections IVII 
and IVII1 we discuss the Pauli-suppression and Bose-enhancement of pair-productions rate in the 
presence of thermal electrons and bosons at a given temperature and chemical potential. Finally, 
in Section IVIIIl discuss the possibility that the pair-production rate can be greatly enhanced by 
the self-focusing phenomenon of laser beam propagating in the plasma of electrons and protons. 

II. SEMICLASSICAL DESCRIPTION OF PAIR PRODUCTION 

The phenomenon of pair production in an external electric field can be understood, in the 
historic first-quantized Dirac picture, as a quantum-mechanical tunnelingprocess of electrons from 
the negative-energy Dirac sea to the positive energy conduction band (la . fll] ]. The electric field 
bends the positive and negative-energy levels of the Hamiltonian, leading to a level-crossing and a 
tunneling of the electrons in the negative-energy band to the positive-energy band. Let the field 
vector E(z) point in the z-direction. In the one-dimensional potential energy 



V(z) = -eAo(z) = ej dz'E{z') (1) 
of an electron of charge — e, the classical positive and negative-energy spectra are 



£±(Pz,P±;z) = ±^{cp z ) 2 + c 2 pl + {m e c 2 ) 2 + V(z), (2) 

where p z is the momentum in the ^-direction, pj_ the momentum orthogonal to it, and p±_ = \p±\- 
For a given energy £, the tunneling takes place from z- to 24. determined by p z = in Eq. ([2]) 

£ = S+(0 3 p ± ;z+) = S-(O lP± ;z-). (3) 

The points z± are the turning points of the classical trajectories crossing from the positive-energy 
band to the negative one at energy £. They satisfy the equations 



V(z±) = TJc 2 pi + m 2 e c* + £. (4) 
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This energy-level crossing £ is shown in Fig. Q] for the Sauter potential V{z) oc tanh(z/£). 




FIG. 1: Positive- and negative-energy spectra £±{z) of Eq. ([2]) in units of m e c 2 , with p z = p± = as a 
function ol z = zjl for the Sauter potential V±(z) (T7T|) for er = 5. This figure is reproduced from Fig. f in 



Rcf. 
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A. WKB transmission probability for Klein-Gordon Field 



In our previous article [12|], using the WKB method we obtain the general formula for the 
electron-positron pair production rate in space- or time-dependent electric fields. In this article, 
using this approach, we further study the electron-positron pair production rate in the presence of 
radiation (photon) and electron fields. First we consider a monochromatic radiation field, which is 
described as an electromagnetic plane wave A^, 

A fi (x) = [A ,A ± (t,z),A g (t)], A = 0, A z = 0, 8^ = ^ = (5) 

with the wave vector k^ = (lj, k), k 2 = 0, propagating in the z-direction k z ^ 0, k± = 0, 

A ± (t,z) = A ± (k)exp^(p2 Z -rt) (6) 

where radiation photon energy-momentum £ 7 = hui, = hk z and = c\pj\. 

The probability of quantum tunneling in the z-direction is most easily studied for a scalar field 
which satisfies the Klein-Gordon equation 

1 2 "1 



ilidp + -A*(z) 
c p 



m y <f>(x) = o, (7) 
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where xq = ct, and the gauge potential A*^(z) is superposition of two parts: the static field 
A^(z) = (Aq(z), 0, 0, 0) of Eq. ([I]) and an electromagnetic plane wave of Eq. ©. Then Eq. (|7|) 
becomes 



ihda + -AJz) 
c 



m y + 2ih-(Ap&>) + e -A 2 4>{x) = o, 



(8) 



where A^A 11 = for the gauge chosen. Since there are only the static electric field in the z-direction 
and electromagnetic plane-wave fields in the perpendicular x, y-directions, both vary only along z, 
we can choose a vector potential with the only nonzero component (fTJ, and make the ansatz: 



= e- i£t / h jv^ h 4> Vx ,e{z)x Px ,e{z), 

a fixed momentum pj_ in the x, y direction and an energy 8, and Eq. ([8]) becomes simply 

d 2 1 e e 2 

{ - ft 2 ^ +P± + m 'c 2 - ^ [S - V(zf + 2^(A ±P± ) - ^A 2 }cb p± , £ (z) Xp± ,£ (z) = 0. 

By expressing the wave function <j> p±t £(z) as an exponential 

^ £ (z)=Ce iS ^l\ 

where C is some normalization constant, we impose the wave equation for (f> p±j g(z) 

1 



(9) 



(10) 



(11) 



<P P± ,s{z) = 0, 



(12) 



which becomes a Riccatti equation for S p 



ihd 2 z S p± , £ (z) + [d z S p± , e (z)] 2 -pl(z) = 0. 



(13) 



where the function p z (z) is the solution of the equation 
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p 2 z (z) = ^[£-V(z)] 2 -p 2 ± -m 2 e c 2 . 



(14) 



The solution of Eq. (|13p can be found iteratively as an expansion in powers of fi: 

(15) 

'lie WKB approximation 



Neglecting the expansion terms after 5^ g (z) = — log p z /z (z) leads to t 
for the wave functions of positive and negative energies can (see e.g. 20 



1/2/ 



2l|) 



i WKB / \ _ 
( P P± ,£ \ z ) - — 



c js™ e (z)/n 



Pz (z) 



(16) 



where £ (z) is the eikonal 



s^Az) = I Pz (z')dz 



p±,£ 



(17) 



Following the differential equations (fTU|) and (fT2|) , the differential equation for the function Xp ± ,e {<fi) 
is given by 

d e 



X P± ,e(z) = 0. 



where we use d^F = k^F', d^F = k 2 F" = 0, k^A^ = 0, A z = and 



Therefore Eq. (fT8|) becomes 



' dz 



-2ihp z — + 2^(A ±P± ) - — Ai 



X Px ,£(*) = o. 



(18) 



(19) 



(20) 



The integral of this equation is 



-Ai-—(A ±P± ) 



dz', 



2c 2 p z cp z 

where the integration constant is absorbed into C in Eq. (jlip . The final WKB-solution is then 



(21) 



■exp+^ / V z (z?)dJ, 



where 



7? 2 (z)=p 2 (z) + 



■(A±p. 



(22) 



(23) 



2c 2 p z cp, 

Between the turning points z^ < z < z+, whose positions are illustrated in Fig. [1] the momen 
turn p z (z) is imaginary and is useful to define the positive function 



1 



k z (z) = dpi + m 2 e c 2 --[£- V(z)] 2 > 0, 



and we define 



K-z(z) = K z (z) 



■(A ±P± ) 



(24) 



(25) 



2c 2 K z (z) CK z (z) 

which depends on the electromagnetic field Aj_. The tunneling wave function in this regime is the 
linear combination 
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Outside the turning points, i.e., for z < z_ and z > z+, there exist negative-energy and positive- 
energy solutions for £ < £ _ and £ > E+ for positive p z . On the left-hand side of z~, the general 
solution is a linear combination of an incoming wave running to the right and outgoing wave 
running to the left: 



exp 



n 



exp 



t / V z dz 



On the right hand of z+, there is only an outgoing wave 

r 



The connection equations can be solved by 



exp 



V z dz 



z + 



C = 0, C± = e ±i7r/4 C/2, T = C+exp 
The incident flux density is 



1 /- 2 + 



K. z dz 



2m P i 



>d z <f> - (d z f) 



Pz /* 

— ( 

me 



which can be written as 



(27) 



(28) 



(29) 



(30) 



(31) 



where f z (z) = p^(^) /m e is the velocity and n—{z) = <f>*{z)<j){z) the density of the incoming particles. 
Note that the z-dependence of v z (z) and ri-(z) cancel each other. By analogy, the outgoing flux 
density is \T\ 2 /m e . 

B. Rate of pair production 



From the above considerations we obtain for the transmission probability 

transmitted flux 



the simple exponential 



Wwkb 



Wwkb(p±,£, A_ 



incident flux 



exp 



n 



z + 



JC z (z')dz' 



(32) 



(33) 



In order to derive from (I32p the total rate of pair production in the electric field we must multiply 
it with the incident particle flux density at the entrance z_ of the tunnel. The particle velocity at 
that point is v z = d£/dp z , where the relation between £ and Z- is given by Eq. (|4|): 

£-V{z-) 



1 



V( C P-L) 2 + m e c4 



(34) 
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This must be multiplied with the particle density which is given by the phase space density 
d 3 p/(27rh) 3 . The incident flux density at the tunnel entrance is therefore 

72 „ 



d£ d 2 p± dp z _ f d£ d 2 p±_ 
dp z (2vr?i) 2 2vr/i ~ s ] M(M)2' 



(35) 



2J 



and the extra factor D s is equal to 2 for electrons with two spin orientations 
It is useful to change the variable of integration from z to (,{z) defined by 

e - v(z) 

and to introduce the notation for the electric field E(p±,£; £) = E[z(p±, £; £)], where z(p±,£; £) is 
the inverse function of (i36l) . the equations in (jlj) reduce to 



(36) 



(37) 



In terms of the variable £, the WKB transmission probability ([33|) can be rewritten as 



Wwkb(p±,£, A_l) = expj 



+ 



-En 



7T^ 

E() 



1 + 



(cpj 



m?c 4 



eA 



;(cp. 



G(p±,f) + 
ff(p±,f)}, 



(38) 



where the photon amplitude A^ is independent of space coordinate. Here we have introduced a 
standard field strength Eq to make the integral in the exponent dimensionless, which we abbreviate 
by 

2 r l 



G( P± ,£) 



nj-i ^E(p ± ,£;()/E 



irJ-i ^E( P± ,£;()/E ' 
The first term in the exponent of (f38|) is equal to 2E c /Eq, where 



E c = mlc 3 /eh 



(39) 
(40) 

(41) 



is the critical field strength which creates a pair over two Compton wavelengths 2Ac = 2h/m e c. 

At the semiclassical level, tunneling takes place only if the potential height is larger than 2m e c 2 
and for energies £ for which there are two real turning points z±. The total tunneling rate is 
obtained by integrating over all incoming momenta and the total area V± = J dxdy of the incoming 
flux. The WKB-rate per area is 



TwKB 



D, 



d£ f d 2 p± 



2TrhJ {2irh)' 



W / wkb(rl,£, A_ 



(42) 
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Using the relation following from ([3 



d£ = eE(z^)dz^, 



we obtain the alternative expression 

Twkb - f dz_ 



d 2 p± 



eE(z-)WwK B \p±,£(z-), Aj.], 



(43) 



(44) 



V± S J 2-KhJ (2ttK) 2 

where £{z~) is obtained by solving the differential equation ()43D . 

The integral over p± cannot be done exactly. At the semiclassical level, this is fortunately 
not necessary. Since E c is proportional to 1/h, the exponential in (|38[) restricts the transverse 
momentum p± to be small of the order of \/h, so that the integral in (|44p may be calculated from 
an expansion of G(p±,£) and H(p±,£) up to the order p\: 

2 ri 



G(p±,S) 



7T 



i ^E(0,£;()/E 



1 



ldE(0,£,()/d( 
2 £(0,£,C) 



c*+... 



G(0,£)+Gj(0,f)5 + ... 



and 



fl-(p±,£) ^ 



1 



ldE(0,£,C)/dC 
2 £(0,5,C) 



where 5 = <5(pj 



= #(0,£) + #5(0, £)£ + ..., 
(cp ± ) 2 /(m 2 c 4 ), 
Gi(0,f) = 



i r 1 c 2 



^T^E(0,£,()/E ' 



and 



H S (0,S) 



/ 1 ^ %m7^/^ o,g;C) 

vrj-i E 2 (0,£;Q/E o 

Into n-lT C 2 < 

2 1 ' J Tri-i (I-C 2 ) 3 / 2 £7(0, CV^o " 



We can now perform the integral over pj_ in 



approximately as follows: 



d 2 p± 



(2vrft) : 
nEc 

Eq 

mlc 2 



exp 



{ 



ttE c 
En 



(l + 6)[G(p,£)+G s (0,£)5\ + 



(i 2 _ a±< 5i/2)[# (0) £) + fr,(0,5)<5]} 



•25&[G(0,£)-iaiH(0,£)] 



4vrft 2 

4^cG(0, 5) 



d(5 e -' r ( E c / £; o)[ 5 G'(o,£:)+<5 1 / 2 a ± H(o,£:)] 



(45) 



(46) 



(47) 



(48) 



(49) 
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where 



a± = e A± ; (50) 
m e c z 

and 

(5(0, 5) = G(0, 5) + G S (Q, £) - \a\H s (Q, 8) (51) 

/ vrEc\ ai£Ml . ( 52) 

The electric fields E(pj_,£; £) at the tunnel entrance z_ in the prefactor of (|44p can be expanded 
similarly to first order in 5. If denotes the solutions of (|34|) at p± = 0, we see that for small 5: 



so that 



Here the extra term proportional to 5 can be neglected in the semiclassical limit since it gives a 
contribution to the prefactor of the order h. Thus we obtain the WKB-rate (|44p of pair production 
per unit area 

Twkb _ I , d z T wkb ( z ) 



V± J V± 

~ D [ d ~ e2EoE ^ c -^[G(0,e)-^H(0,£)] 
S J 87r s h 2 cG(0,£{z)) 

• {l + Tr^e^tl + Erf^)]}, (55) 

where z is short for z®_. At this point it is useful to return from the integral / dz-eE(z-) introduced 
in (|44p to the original energy integral / d£ in (|42p . so that the final result is 
Twkb _ f ^wkb(^) 



Vj_ J V ± 

~ D eE ° f d£ . 1 c -ffi[G(Q.g)-Mg(o.g)] 



4vr 2 nc7 27rftG(0,£) 
• {l + 7T 1/2 ?9e 1?2 [l + Erf (0)]}, (56) 

where ^-integration is over all crossing energy- levels. 

These formula can be approximately applied to the 3-dimensional case of electric fields E(cc, y, z) 
and potentials V(x, y, z) at the points (x, y, z) where the tunneling length a = z+ — z_ is much 
smaller than the variation lengths 5x±_ of electric potentials V(x, y, z) in the xy-plane, 

1 1 sv 

- > - — . 57 

a V ox±_ 
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At these points (x,y,z), we can arrange the tunneling path dz and momentum p z (x,y,z) in the 
direction of electric field, corresponding perpendicular area d 2 V± = dxdy for incident flux and 
perpendicular momentum p^. It is then approximately reduced to a one-dimensional problem in 
the region of size 0(d) around these points. The surfaces z_(x_,y_,<5) and z + = (x + ,y + ,£) 
assciated with the classical turning points are determined by Eqs. (|36h and Eqs. (I37|) for a given 
energy £. The WKB-rate of pair production (I55p can then be expressed as an volume integral over 
the rate density per volume element 

Twkb = / dxdydz d ^ ^ B = / dtdxdydz f . (58) 

J dx dy dz J dt dx dy dz 

On the right-hand side we have found it useful to rewrite the rate Twkb as the time derivative of 
the number of pair creation events dNyjKB/dt, so that we obtain an event density in four-space 



ri 4 iVWKB D e 2 E E(z) -lE±lG(0,E)-\a\H(S)£)] 

dtdxdydz ~ s 87r 3 hG(0,£(z)) e 

• {l + 7r 1/ W 2 [l + Erf (■&)]}, (59) 

Here x,y and z are related by the function z = Z-(x,y,£) which is obtained by solving ([43 p . 

It is now useful to observe that the left-hand side of (|59j) is a Lorentz-invariant quantity. In 
addition, it is symmetric under the exchange of time and z, and this symmetry will be exploited 
in the next section to relate pair production processes in a z-dependent electric field E(z) to those 
in a time-dependent field E{t). 

Attempts to go beyond the WKB results (|55|) or ()56[) require a great amount of work. Corrections 
will come from three sources: 

I from the higher terms of order in {fi) n with n > 1 in the the expansion (|15p solving the 
Riccati equation ()13|) . 

II from the higher terms of the perturbative evaluation of the integral over in Eqs. ()42|) or 
(I44p when going beyond the Gaussian approximation. 

Ill from perturbative corrections to the Gaussian energy integral (|56p or the corresponding 
z-integral ([53]) . 

All these corrections contribute terms of higher order in h. 
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C. Including a Smoothly Varying B(z)-Field Parallel to E(z) 

The above results can easily be extended to allow for the presence of a constant magnetic field 
B parallel to E(z). Then the wave function factorizes into a Landau state and a spinor function 
first calculated by Sauter [l|. In the WKB approximation, the energy spectrum is still given by 
Eq. (J2J), but the squared transverse momenta p\ is quantized and must be replaced by discrete 
values corresponding to the Landau energy levels. From the known nonrelativistic levels for the 
Hamiltonian p 2 L /2m e we extract immediately the replacements 



c 2 p\ = 2m e c 2 x ( J — > 2m e c 2 x 



n = 0,1,2,- (60) 



where g = 2 + a/ir + . . . is the anomalous magnetic moment of the electron, ujl = eB/m e c the 
Landau frequency, with a = ±1/2 for spin- 1/2 and a = for spin-0, which are eigenvalues of the 
Pauli matrix a z . The quantum number n characterizing the Landau levels counts the levels of 
the harmonic oscillations in the plane orthogonal to the z-direction. Apart from the replacement 
(|60p . the WKB calculations remain the same. Thus we must only replace the integration over 
the transverse momenta / d 2 p±/(2-rr?i) 2 in Eq. (|49[) by the sum over all Landau levels with the 
degeneracy eB/(2irhc). Thus, the right-hand side becomes 

eB ~7r(Ec/g )[G(0,g)-f a j g(0,g)] -7r(iJ/i?o)[(n+l/2+g ( T)G(0,£:)+(n+l/2+9 ( 7) 1 / 2 a ±J H'(0,£:)] 

The approximate result is, for spin-0 and spin-1/2: 

a 2 /^n r , e-^/^)[ g (°- £ )-^l g ( - £ )l/o, 1/2 (gG(0,g)/£b) (62) 
47r^n,cG(0, t) 

where 



/0(x) = sinh^' h / 2{x) = 2 ^nh^ C ° Sh — (63) 



In the limit B -)• 0, Eq. ([63]) reduces to Eq. (l49j) . 

The result remains approximately valid if the magnetic field has a smooth z-dependence varying 
little over a Compton wavelength A e . 

In the following we shall focus on nonuniform electric fields without a magnetic field. 



III. TIME-DEPENDENT ELECTRIC FIELDS 



The above semiclassical considerations can be applied with little change to the different physical 
situation in which the electric field along the z-direction depends only on time rather than z. Instead 
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of the time t itself we shall prefer working with the zeroth length coordinate xq = ct, as usual in 
relativistic calculations. As an intermediate step consider for a moment a vector potential 



A^ = (A {z), 0,0, A z (x )), 



(64) 



with the electric field 



E = -d z A (z) - d A z (x Q ), x = ct. 



(65) 



The associated Klein-Gordon equation ([7]) reads 



ihd + -Aq(z) 



+ h z d: 



2fj2 



iW z + -A z (x Q ] 



m 2 e c 2 + 2ih-(A fl d' 1 ) + -^Ai } <j>(x) = 0.(66) 



The previous discussion was valid under the assumption A z [xq) = 0, in which case the ansatz 
led to the field equation (llOp . For the present discussion it is useful to write the ansatz as 



[x = e 



-ip x /h ip ± x±/h 



^Px,Po( 2 )X P± ,po(^) 



(67) 



with po = £/c, and Eq. (|10p in the form 



8 -el dz 1 E(z') 



1 2 



d 2 e e 2 I 

pi - rale 2 + h 2 ^ - 2-(A_lp_l) + ^A 2 ± <f> P± , Po (z)x Px , Po (z) = ($68) 



Now we assume the electric field to depend only on xq = ct. Then the ansatz: 



( x \ — e ipzz/h e ip ± x ± /h 



(69) 



leads to the field equation 
( d 2 



2„2 



dx< 



2~P±~ m e C 



dx' E(x' ) 



2-(A_lp_l) + — A^ \ 4> P± ,p z (x )x P± ,pA x o) 



If we compare Eq. (|70|) with (|68p we realize that one arises from the other by interchanging 



zHiq, p± — > ip±, c — )• ic, E — > —iE. 



(71) 



With these exchanges we may easily calculate the decay rate of the vacuum caused by a time- 
dependent electric field E{x$) using the above-derived formulas. 



IV. APPLICATIONS 

In this section, we apply formulas (156|) or (I55p to two typical external field configurations 
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A. Step-like electric field 

First we check our result for the original case of a constant electric field E(z) = cEq where the 

potential energy is the linear function V(z) = —eE§z. Here the functions (I39p and (140p become 
trivial 

G(0,£) = - f d(Jl-e = l, G s (0,£)=0, (72) 
vr J -i v 

H(0,£) = - [ l d((Jl - er 1 = 1, H S (0,£) = -1, (73) 

(5(0,5) = 1 + 4/2 (74) 

which is independent of £ (or z_). The WKB-rate for pair-production per unit time and volume 
is found from Eq. (|55|) to be 

% - ^ 8 Sii + i i/2 e "^ (1 " ai/2) + ^ i/2 ^ 2 [i + Erf( ^ )] } • (?6) 

where V = dz_ V± . We find that photon field amplitude squared a\ ([50]) gives rise to an exponential 
factor of enhancement e^ nEc ^ E °^ a2 L^ 2 \ we will turn to this point later. For a±_ — > and •& — > 0, 
Eq. ([76]) goes to the correct expression found by Sauter 1|, Heisenberg and Euler Q|, and by 
Schwinger 

B. Sauter electric field 

Let us now consider the nontrivial Sauter electric field concentrated to a thin slab in the xy- 
plane with a width I in the ^-direction. A field of this type can be produced, e.g., between two 
opposite charged conducting plates. The electric field E(z)z in the z-direction and the associated 
potential energy V(z) are given by 

E(z) = £ /cosh 2 (z/£) , V{z) = -a m e c 2 tanh (z/£) , (77) 

where 

a = eE £/m e c 2 = (£/ X C )(E / E c ). (78) 

From now on we shall use natural units in which energies are measured in units of m e c 2 . Figure [1] 
shows the positive and negative-energy spectra £±{z) of Eq. ([2]) for p z = p± = to show the 
energy-gap and energy-level crossings. From Eq. ([4]) we find the classical turning points 

z± = l arctanh^^ 7 = '- b.'^^. (79) 
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Tunneling is possible for all energies satisfying 



VTTS + a>£> y/T+S - a, 



(80) 



for the strength parameter a > y/1 + 6. 

We may invert Eq. (|36p to find the relation between £ and z: 



< e r\ ♦ h £ + CVT±l £, a + £ + (VTT6 
z = z(p± : c; C) = t arctann = - in . 

In terms of the function z(p±,£; £), the equation (j79|) . reads simply z± = z(p±,£; ±1) 
In Ref. jl^ . inserting (|8ip into the equation for E(z) in Eq. ()77|) . we obtain 
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(81) 



= E 



1 



= E(p ± ,£;(), 



(82) 



and calculate G(0,£) and G s (0,£) of Eqs. ([39D, $M and (S3) : 

r i 1/2 r 

G(0,£) = 2a 2 - a {a - £) 2 - 1 - a (a + £) 2 - 1 



(83) 



and 



G(0,£)+G*(0,£) 



-1/2 



+ 



(a + £) 2 -l 



-1/2 



(84) 



The integral over all energy-level crossings permitted by the energy inequality (|80p is dominated 
by the region around £ ~ 0, where the tunneling length is shortest [see Fig. Q] and tunneling 
probability is largest. Both functions G(0,£) and Gs(0,£) have a symmetric peak at £ = 0. 
Around the peak they can be expanded in powers of £ as 



G(0,£) = 2[a 2 - a{a 2 - 1) 1/2 ] + 



cr 



(a 2 - l) 3 / 2 



£ 2 + 0{£ 4 ) 



G (a) + -G 2 (a)£ 2 + O(£ 4 ), 



(85) 



and 



G(0,£) + G s (0,£) 



a 



( a 2- 1)1/2 2 ((J 2 -l) 5 / 2 

G (a) + ^G 2 ((r)f 2 + C»(f 4 ) 



(86) 



The exponential e - wG (°> £ ) E c/ E o k ag a Q auss i an peak around £ = whose width is of the order of 
1/E C , which in physical units is a first-order term in h. This implies that in the semiclassical limit, 
we may perform only a Gaussian integral and neglect the ^-dependence of the prefactor in (|56|) . 
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Analogously, the function H(0,£) of Eq. (|46p is given by 

fi 



H{0,£) 



1 i 



2a 



+ 



a(l - 2a' 



(1 + (72)1/2 (l + ff 2)5/2 
Io(ff)+H 2 ((T)f 2 +0^ 4 



£ 2 + 0(£ 4 



(87) 



and H s (0,£) Eq. (08 



^(0,5) = ~H(0,£)-- [ d(— 

2 7T J-l (1 



C 2)3/2 1 _/ i _ £ y 



3 r 



The equations ([51]) and (152]) become 



H (a) + H 2 (a)£ 2 +0(£ A ). 



1^, 



(88) 



G(0,£) = G (<x) + ^G 2 {a)£ 2 



+ ^1 



tf 2 = 



H {a) + H 2 {a)£ 2 ] + 0{£ A ) 
ttE c \ a 2 ± H 2 (0,£) 



(89) 



E J 4 G(0,f) ' 



So / 4 G (a) + f ai^o(^; 



+ o(£ 2 ). 



(90) 



Recalling that £ in this section is in natural units with m e c 2 = 1, we must replace / d£ in 
the pair-production rate (|56p by m e c 2 j d£ and can perform the integral over £ approximately as 
follows 

eEom e c 2 1 



Twkb 
V± 



4vr 2 nc Go + %a 2 H 



e -7v(E c /E )(Go-^alH ) j d£ r -7z(E c /E )(G2-a 2 ± H2) S 2 /2 



D R 



eE 



2-kU 



■ir{E c /E ){G -\a\H ) 



(91) 



'47r 2 7ic G + \a\Ho 27rft[(G 2 - a 2 ± H 2 )E c /2E } 1 / 2 ' 
For convenience, we have extended the limits of integration over £ from the interval (— 1 + cr, 1 — a 
to (—00,00). This introduces exponentially small errors and can be ignored. 



COHERENT AND THERMAL STATES OF PHOTON FIELDS 



Now we turn to the calculation of the average of = (eA ± ) 2 /(m 2 c i ) [Eq. ([50])] of photon fields 



in coherent monochromatic state and thermal state. 
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The pair-production rate (|59p depends on the transverse amplitude A±(t,z) ([5]) of monochro- 
matic electromagnetic wave cu = |k| = k z , we need to take the average 

rf 4 iVwKB 
dt dx dy dz 



over amplitudes Aj_. Using the convexity inequality |23J] 



(e A I)>e^, (92) 



we can obtain the low bound in the pair-production rate Eq. (1591) , 

For the case of monochromatic electromagnetic wave © with its transversed amplitude A±(k), 
= k), and the corresponding electric component 

E±(<, «) = -qI A ^^ z ) = -*^Aj.(<, z), (93) 

and the maximal amplitude Speak = w l Aj_(&)|/c. In a monochromatical state of laser photon, 
averaging over one period T = 2it/oj, we have 

(Ai(M)> = ^ [ T dtAl(t,z)A ± (t,z) = A* x (k)A ± (k), (94) 

and 

1 2 l/ (eA ± ) 2 \ 1 e 2 m 1 /m e c 2 \ 2 /Spe^X 2 

2 ( ^ > = 2\^^/ = 2^ A ^ )A±(fc) = 2l v ^rJ l^rJ • (95) 

This shows the range of -Ep ea k and fro; in order to have a significant enhancement. 
Let us now consider a general gauge field (h = c = 1), 

A±(x) = y ^M^A^AOe^ = / j^^6 + (k 2 )A ± (k)e ik *, (96) 
where is the dispersion relation of electromagnetic fields A±(x). Averaging over space-time, we 



have 



(Al(x)A ± (x)> = J d*xA* ± (x)A ± (x) = J Al(fc)A ± (fc) 

d 3 k 

-A* ± (k)A ± (k), (97) 



where ojk = \k\. Here A* L (k)A±(k) is the number density of photons in the momentum state k. In 
the case of thermal photons at temperature T, it is equal to 
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Here and in the sequel T will be measured in natural units in which the Boltzmann constant Ub is 
equal to unity. Hence the finite-T version of Eq. (I95p is 



1, 3 2a f d 3 k 

2 M = ^JWN~k m 



2a f d 3 k 1 a ( T 



2m2 c 4 J (27r) 3 a; fc e^/ T - 1 12 \m e c 2 
This shows that the enhancement is very small for T < m e c?. 



(99) 



VI. THE PRESENCE OF ELECTRONS WITH TEMPERATURE AND CHEMICAL 

POTENTIAL 

In the previous section (|V|) . we have considered the pair-production rate in the case of thermal 
photons at temperature T (|98p . We further consider the pair-production rate in the gas of thermal 
electrons at the same temperature T. These thermal electrons are in the Fermi distribution, 

fe(£e,Ve,T) = e{£e _^ /T + i , (100) 

where the electron energy-level 

£e = [{cp e f + mlc 2 } 1 ' 2 , (101) 

the associated electron number-density is 

n e (m e , Me ,T) = 2 J ^ ^— ^ — , (102) 

and the chemical potential [i e > that is related to the total number M e of electrons. The rate of 
pair-production is given by 

Twkb [ d£ f d 2 p A 

where energy-level-crossings for pair-productions, and the inserted Pauli-blocking factor [1 — 
f e (£ , fj, e , T)]s = £ e gives a probability whether the energy-level £ e = £ is occupied. This limits the 
phase-space permitted by available energy-level-crossings "£" for pair-productions. 

In the case of constant electric fields (|76[) . the pair-production probability Wwkb(pj_, £, A) 
is independent of the energy-crossing- level u £" . We first consider the constant electric field Eq 
confined within the finite "box" region [— £/2,£/2] in the z-direction, and the range of energy-level 
crossing is and £_ > £ + , where 

£± = TeEol/2, (104) 



D.J ^[l-/ e (^/ie,T)] J J^WwsB(px,e,A), (103) 
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£_ > and £j r < (Fig.[T]presents a similar case). Furthermore, we assume that in this finite "box" 
region there are electrons whose density is given by n e (fJ, e ,T) [see Eq. (I102p ], We can calculate the 
pair-production rate by integrating over energy-crossing-levels 

f £ - d£ T ( e £ -l T + e^l T \ 

k m*-M<^ = m*\ wX+J ) < 105) 

where £_ > \i e > £+. As a result, the pair-production rate per area (j76[) is modified as follows 

Elk n ^L^ ( e £ -l T + e^ T \ 
V ± ~ s 8^h 2 c \ e £ +/ T + e^/ T J 

x ' 2W0 e-^ (Mal)/2) (l + vrVV 2 [l + Erf(*)]} , (106) 

where i? = i9((a^_)) [see Eq. (f75]) ]. Eq. (fT05|) plays a suppression factor in the pair-production rate 
(fT06l) . 

In the low-temperature limit T j /i e <C 1 and electron chemical potential goes to its Fermi energy 
Ep, /i e — > Ep the leading order of Eq. (|105p is given by 

g--/j e ^ eE £/2 - Ep 
2-nfi 2ttH 1 ' 

When E- — ix e = eEo£/2 — Ep = 0, which means all energy crossing-levels for pair-productions 
are fully filled by electrons, it leads to a complete Pauli-blocking, and vanishing of pair-production 
rate (I106j) . In the high-temperature limit T j /x e 3> 1, the leading order of Eq. fll05f) is given by 

and Eq. ()lQ6j) correctly goes back to the expression (fT6j) . 

In the case of Sauter electric field (|77p for the semiclassical limit, neglecting £ -dependence in 
the prefactor, we consider the following Gaussian ^-integration (see Fig. [1]) 

where the exponential factor plays a cutoff at 

£± = ± i^j' 2 (£ c /£o)- 1/2 (G 2 - alH 2 r 1/2 , (110) 

a > £+ and £_ > — a. To see the Pauli-blocking effect, we further neglect ^-dependence of the 
exponential factor of Eq. (|109p . and assume the maximal pair-production probability at £ = 0, 
and Eq. (|109p approximately becomes 

a-l d £ T f p{a—l)/T I n e /T\ 
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In consequence, the pair-production rate in the Sauter field Eq. (|9Tj) multiplied by this expression 
factor (|llip . We have the same discussions on the high- and low- temperature limits by the 
replacement a — 1 =4> eE£/2. 

VII. THE PRESENCE OF BOSONS WITH TEMPERATURE AND CHEMICAL 

POTENTIAL 

In the previous section (IVip . we have considered the suppression of pair-production rate in the 
presence of thermal electrons at temperature T and chemical potential [i e . Considering charged 
bosons 4>, we further consider the enhancement of pair-production rate in the gas of thermal bosons 
at the same temperature T. These bosons are in the Bose-Einstein distribution 

U(S 4t ,fi <f „T) = e( ^_^ )/T _ 1 , (112) 

where the boson energy-level 

^ = [(^) 2 +m^] 1/2 , (113) 

the associated boson number-density 

n^m^T) = 2 J ±2jL ^_* )/T — , (114) 

where £ ( f > > fi^ and the chemical potential fi^ is related to the total number N$ of bosons. The 
rate of pair-production is given by 

Tj VT = ° S I ^[ l + U(^H^)] J -0^W WK b(p±,£,A), (115) 

where "£" energy-level-crossings for pair-productions, and the inserted Bose-Einstein enhancement 
factor [1 + f,p(£ , fx,p,T)]g = £, gives a probability that more particles can occupy the energy-level 
£fj, = £. This enlarges the phase-space permitted by available energy-level-crossings "£" for pair- 
productions. 

In the case of constant electric fields (|76p . the pair-production probability Wwkb(rl> -4) 
is independent of the energy-crossing-level We first consider the constant electric field Eq 

confined within the finite "box" region [— £/2,£/2] in the 5-direction, and the range of energy-level 
crossing is [£+,£_] and £L > £+, see Eq. (|104p . Furthermore, we assume that in this finite "box" 
region there are bosons whose density is given by (|114p . We can calculate the pair-production 
rate by integrating over energy-crossing-levels 

r £ - d£ r , . ,„ „„ 1 ,„ _ T . (e £ ~/ T -e^ 
>£+ 



r S- d£ 1 T ( e £ -/ T - e^/ T \ 
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where (£L — 8 + ) = eE$l and /i<£. As a result, the pair-production rate per area ([76]) is 

modified as follows 



r WKB _ D eE o T 



V± " 8tt 3 ^ 



2, 



C 



2 / ^T^VT N 
-(g_-g + )-ln( £+/T _^ /T 



x - / 2 . e -^f(M4)/2) (! + V 2 [! + ErfW] \ . (117) 
1 + (a_|_)/2 

Eq. ()116p plays an enhancement factor in the pair-production rate f)117[) . 

In the low-temperature limit T/j/x^j <C 1, the distribution (|112|) shows bosons undergo the 
Bose-Einstein condensation, by going to the energy level 8^ = fi^, and the momentum states 
P(j> = ~ rn <t>c l > the leading order of enhancement factor in Eq. (|1 17j) is given by 

^(S- - 8+) = ^(eE t) > 0. H . (118) 

We find that the enhancement factor is two, by comparing Eq. (j!17|) with Eq. (|76|) . In the high- 
temperature limit T 3> 1 and T < {8- — 8 + ) the leading order of Eq. (|116p is 

h8^-8 + )-\n £ -^± = l( e E £)-ln^^, (119) 

and we find that the enhancement factor is 

K 1 -^ 1 "^)' (120) 

by comparing Eq. (|117p with Eq. (fT6j) . 

In the case of Sauter electric field (j77|) for the semiclassical limit, neglecting ^-dependence in 
the prefactor, we consider the following Gaussian ^-integration (see Fig. [T]) 

Cm^'^ J '^ + >•*• r » - < £ - - f +> - in - ^ 

where the exponential factor plays a cutoff given by Eq. (jllOj) . To see the Bose-Einstein enhance- 
ment, we further neglect ^-dependence of the exponential factor of Eq. (|12ip . and assume the 
maximal pair-production probability at 8 = 0, and Eq. (1121 1) approximately becomes 
ra-\ A£ T ( {*-\)lT _ H IT\ 



In consequence, the pair-production rate in the Sauter field Eq. (|91|) multiplied by this enhancement 
factor (|122p . We have the same discussions on the high- and low- temperature limits by the 
replacement a — 1 =4> eE£/2. 



22 

VIII. THE PRESENCE OF A NEUTRAL PLASMA OF ELECTRONS AND PROTONS 

In the previous section, we have considered the pair-production rate in the presence of thermal 
photons at temperature T (198|) . Another physically interesting environment is the presence of 
a neutral plasma composed of electrons and protons. The two charge components can oscillate 
against each other and modify the electric field available for pair creation. For simplicity let us 
assume the protons to form a charged lattice and let us ignore the temperature Toebyc associated 
with the lattice phonons. The electrons are distributed in the lattice so as to screen electric fields of 
the proton charges and of the external electric potential Aq. In such an equilibrium configuration, 
we shall assume the electrons to be in a thermal equilibrium at a temperature T and chemical 
potential (M e , so that their Fermi distribution is by Eqs. (jlOOp and (jlOip . The associated electron 
number-density (|102p . energy-density and pressure 



The chemical potential fi e > is fixed by the total number Vn e (m e , /j, e ,T). The gradient of 
electron-gas pressure balances all electric forces. 

Due to perturbations, these electrons deviate from their equilibrium positions, and this may 
lead to the coherent plasma oscillation of electrons in the proton lattice. In order to study this, we 
first neglect the dissipative terms, and describe perturbation of these electrons as a simple perfect 
fluid, whose energy-momentum tensor, 



where the flat metric g^ u = (—,+,+,+) and electron four velocity U£. In the energy-momentum 
tensor (|125p . 5n e , 5e e and 5P e are perturbations of proper number, energy densities and pressure 
in comoving frame of electron fluid. Such plasma oscillation of electrons around the equilibrium 
configuration in the proton lattice can be described by the continuity equation, energy-momentum 
conservation, and the Maxwell equations yield 




(123) 



(124) 



STf = SP e g^ + (SP e + 5e e )U?U- 



rV 

e ' 



(125) 




(126) 




(127) 




(128) 
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where 5F^ y is the strength of fluctuation electromagnetic fields due to the fluctuating electric four 
current 



6J>* = e(5n p U^ - 5n e U^). 



(129) 



Herera p is the proton number-density, and Up = (1,0,0,0) four-velocity of the protons On the 
r.h.s. of Eq. (|127|) . the dissipative term 



eSnpU^SF^, 



(130) 



indicates an Ohmic heating 5Q, and we assume that this term is negligible for 5n p ~ 0, SQ = 
5S/T ~ and the entropy S is approximately conserved (5S ~ 0). This is consistent with non 
dissipative energy-momentum tensor (|125p we have adopted for electrons. In consequence, the 
energy-momentum conservation along four-velocity U£, i.e., U{?d u (5T e u ^) = 0, gives the first law of 
thermodynamics in the form 



5P e + 5e e = [i e bn e 



(131) 



corresponding to the equation of state 8P e /5e e = k? = const, for a isothermal process of constant 
temperature T. 

As discussed, the electrons deviate from their equilibrium positions, thereby creating a small 
electric potential 5Aq, associated with the fluctuating electromagnetic field 5F^ in Eq. (|128p . The 
perturbed electron distribution / e (£ e , /i e , T, 5Aq) is given by the replacement 



£ e —?■ £ e - e5A , 



(132) 



in the electron distribution (jlOOp . Expanding perturbed electron distribution / e (£ e > /^e, T, 5Aq) up 
to the leading order 5Aq, we obtain 



f e (£ e ,H e ,T, 5A ) w f e (S e , n e , T) 
and an electron number-density fluctuation 
6n e (m e ,fi e ,T) 



l + ^5A e( £ °-^ T f e (£ e ,Ve,T) 



(133) 



,(£ e -fle)/T 



f d3 P* 
T / 



:SAn 



T 



as well as energy-density fluctuation 



(2vrfr) 3 [ e (£ e -Me)/T + 1 y 



n e (m e ,fi e ,T) - 2 J —j- 



(2tt) 3 [ e (£e-ne)/T + !]2 



(134) 



2e f d 3 p e £ e e^ £e ~^/ T 

5e e (m e , u e ,T) « — 5Aq / - — — ^ r> — e5Aon e (m e , u e , T) 

y ' ; T J (2vr/i) 3 [ e (£e-M/T + i] 2 K ' 



T 



z5An 



e e (m e , fi e , T) - Tn e (m e , fJ, e ,T) - 2 J 



d 3 p e 



(2vr) 3 [ e (£e-Me)/T + !]2 



.(135) 
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This yields an electron pressure fluctuation 

5P e (m e ,n e ,T) « 2e5A 



d 3 p e 



-{£ B -He)/T 



(2tt^) 3 [ e -(^-Me)/T + !] 
e5A n e (m e , He,T), 



(136) 



which propagates through the electron gas. 

In order to study the propagation of such plasm, we consider the Maxwell equation (I128j) for 
the fluctuation field 5Aq caused by the charge fluctuations 



V 2 5A - 4 ^o^o 



-4ire[5n p — 6n e ], 



(137) 



where the velocity is given in units of the speed of light in vacuum c: 



6Pe 

8e P , 



Tn e (m e ,fj, e ,T) 



d 3 p e 



(138) 



£e{ m e, [J-e,T) Tn e (m e , /i e ,T) 2 J [ e (£ e -n e )/T +1 ]2 

This is a constant k 2 in an isothermal process of constant temperature T. We shall ignore the 
much smaller fluctuations of the proton distribution 5n p ~ 0. Inserting Eq. (|134|) and a plane wave 
ansatz <5Aq = e - lw *+ lkx 5 we obtain an the energy-spectrum for the plasma waves 



^^(IklHo^ + ^kl 2 , 



(139) 



where 



2 _ 



a 2 c 3 



,(£e-fJ. e )/T 



2ira c 



2^3 



7T 2 [ e (£e-Me)/T _|_ 1] 

n e (m e ,fj, e ,T) - 



2 ' 



1 



(140) 



(2vr) 3 [ e (^-fie)/T + !]2 

is the plasma frequency of the electron gas in the proton lattice. These plasma oscillations prop- 
agate through the plasma with a transverse electromagnetic wave A±(x) with two transverse 
polarizations. Their propagator is given by 



dij /cj/cj/|k| 2 



(141) 



w fc" w pi( k ) 

that we call plasmaon field whose energy dispersion is given by (|139|) . corresponding to massive 
photons. Their excitation energies will be in thermal equipartition with the thermal state of 
electrons in the same temperature T. In consequence, the thermal distribution function of these 
massive photons is given by Eq. (|98p with the energy dispersion-relation w 2 |(|k|) of (|139p . Following 
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the same calculations from Eqs. (|97ti99|) . we calculate the average of a\ = (eAj_) 2 /(m 2 c 4 ) [Eq. ([50 
of massive photon fields in thermal plasma state, 

1 , 2 i a f d?k 



2^ ~ J (^p /plW 



a 



(142) 



" 2mic'J (2 I r)3^ 1 (|k|) e "=,(|k|)/r_ 1 - 

For the case that temperature T is much larger than the plasma frequency cj p i, Eq. (|142[) is 
approximately equal to Eq. (|99p . while for the case that T is much smaller than the plasma 
frequency cj p i, Eq. (|142l) is approximately proportional to a(co p ih/m e c 2 ) 2 . 

It is interesting to discuss the case that a monochromatic electromagnetic wave ([6]), (|93p - (195p 
propagates through the plasma of electrons in the proton lattice. Define a dielectric constant 
e = 1 + Xe, where the susceptibility Xe is given by Eqs. (|137|) and (|139|) 

Xe = ~ 2 \ U2 , (143) 
ui z — ir |k| z 

as a function of the frequency oo and wave- vector |k| of the monochromatic electromagnetic wave 
(laser beam) propagating in the plasma. The displacement field strength in the plasma D = eE. 
For large frequencies w 2 /|k| 2 3> v 2 , Xe ~ ~ UJ p\/ uj2 ^ the dielectric constant e w 1 — uj^/uj 2 and e » 1 
for uj 2 3> oj 2 y While for small frequencies w 2 /|k| 2 <C v 2 , Xe ~ + w pi/|k| 2 f 2 , the dielectric constant 
e ~ 1 + a; 2 1 /(|k| 2 v 2 ). The resonance appears at uj 2 = w 2 es = |k| 2 f 2 , at which the dielectric constant 
|e| 3> 1, and the displacement field D greatly increases. 

There are no imaginary damping terms in the denominators of Eqs. (|14ip and (|143p . because 
we use the perfect fluid stress tensor (|125|) for the electron plasma. In particular we neglect Ohmic 
heating in Eq. (|127p . If take this into account in ()130p . we have the following energy dissipation 
per electron in a period T = 2tt/lo: 

<** ~ ~* (£) (#) (m) " " eT (£) ■> (£) " ^ &) (144) 

where 7 ~ 1 is a Lorentz factor. This small dissipative term <5<?diss should be added into Eq. (|132|) . 
namely, replacing the energy perturbation c5Aq by 



e5A + <5£ diss = e5A c 



i+f(27r) (^y 



(145) 



This creates an imaginary damping term in the denominators of Eqs. (|14ip and (|143p . limiting the 
life time of plasmons via a finite width of the resonance. 

However, a great increase of displacement field D at the resonance for co 2 = co 2 es does not yet 
enhance the pair-production rate. The expectation |( a j_) i n Eq. ([95]) for doing this is purely due 
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to electric field E of laser beams, and ultrahigh-intensity laser beams are required. Help can come 
from the self-focusing phenomenon of ultrahigh-intensity laser beams propagating in the plasma 
of electrons and protons. These can be used in principle realize also a large electric field, and thus 
a large term |{tij_) ([95]) . If laser intensities are larger than a certain threshold critical power (see 



review 



m) 



P cr = ^¥-17f^V G W, (146) 

for relativistic self-focusing, the laser pulse can be self-focused when propagating though a plasma 
of electrons and protons with the plasma frequency lo p \ ()140p . It will be interesting to measure 
the electron-positron pair production by a self-focused ultrahigh-intensity laser beam in such an 
environment. 

To end this section we note that if electrons were bosons, one can do calculations by using the 
Bose-Einstein distribution instead of Fermi one. The discussions and conclusions are the same. 
The total pair-production rate receives a factor of suppression and enhancement that are discussed 
in the previous Sees. [VI] and IVHI 



IX. SUMMARY AND REMARKS 

, we studied the process of electron-positron pair production from the vacuum as 
a quantum tunneling phenomenon, we derived in semiclassical approximation the general rate 
formulas ()55[) and (|56p with ^a^_ = 0. They consist of a Boltzmann-like tunneling exponential, 
and a pre-exponential factor, and are applicable to any system where the field strength points 
mainly in one direction and varies only along this direction. In this article, we generalize these 
formulas to the presence of a monochromatic electromagnetic wave ([6]), ()93p - (|95p in addition to a 
static electric field in one direction. We have also considered the system of electrons and charged 
bosons at finite temperature and chemical potential. In several cases, we calculate and discuss the 
factor (^cijj for enhancing pair-production rate. In particular, we consider the plasma of electrons 
and protons, and point out the self-focusing phenomenon of ultrahigh intensity laser beams in the 
plasma possibly gives rise to a larger factor (^aj_) for enhancing pair-production rate, and this 
could be experimentally relevant for observing pair production in laboratories. 

Apart from its purely theoretic interest, these formulas for the pair-production rate in one- 
direction nonuniform fields and laser fields can be relevant for studying the collisions of laser beams 
fl^ and heavy ions 25-27], an explanation of the powerful Gamma Ray Bursts in astrophysics 



5. 
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and strong electric fields in the surface shell of compact star s |29l| . It is also interesting to 



use these formulas to study the phenomenon of plasma oscillations 
after their creation in electric fields. 



301 ] of electrons and positrons 
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